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1. Description of Generator 

Black holes are produced in proton-proton collisions and decay to two- 
particle final states [l| . Electric charge, QCD colour, and spin are conserved, 
while baryon and lepton number need not be conserved. PYTHIA 8 is used 
for the initial- and final-state radiation, multiple interactions, final state par- 
ton showers, fragmentation, and decay processes The QBH generator 
interfaces to PYTHIA 8 using the Les Houches accord (LHA) jij. 

The first black hole generators were TRUENOIR written by Dimopoulos 
and Landsberg^ and an unnamed generator used in the first black hole 
study specific to ATLAS by Tanaka et al. [6]. These generators were made 
quick and efficient for specific studies by averaging over some dynamical quan- 
tities, rather than generating them according to probabilistic distributions. 
The first two true Monte Carlo event generators for black hole production 
and decay were CHARYBDIS [7] and CATFISH jgj]. A comparison between 
CHARYBDIS and CATFISH can be found in Ref. Q- 



The most recent generators CATFISH [8|, Charybdis 2 10], and Black- 



Max all include the facility to generate black holes decaying to two- 
particle final states. It is possible to generate two-particle final states with 
CATFISH but there are some caveats. One sets the final A^-body decay to 
two particles and raises the black hole mass threshold at which Hawking evap- 
oration ends Qmin to be above the minium black hole mass and Planck scale. 



1 lit tp : / /hep . brown, edu /users/ Greg/TrueNoir / 
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If Qmm is set too low, one will occasionally produce black holes above Q m \ n 
resulting in more than two particles in the final state. If Q m i n is set too high, 
one will occasionally produce black holes with high mass that only decay to 
two-particle final states rather than the anticipated Hawking evaporation. 



BlackMax implements the suggestions of Meade and Randall verbatim [12 
The decay species are chosen according to thermal greybody- modified dis- 
tributions, which is perhaps not applicable to black holes near the Planck 
scale Charybdis 2 allows the user to set parameters so that effectively a 
black hole decaying to two-particle final states is generated. CATFISH and 
Charybdis 2 do not take colour explicitly into account when determining the 
two-particle final state. In all three cases, the generators do not allow for 
the ability to generate a particular black hole state, which is a representa- 
tion of SU(2) colour, charge, and spin. They thus act as inclusive black 
hole generators and sum over all possible partons giving rise to black holes. 
In addition, these generators do not allow the ability to easily generate jet- 
lepton or jet-boson topologies. Although baryon number can be violation by 
these generators when run in standalone mode, the partons in the output 
record can not be hadronized by programs that are typically used to provide 
the hardonization, fragmentation, and decay process. QBH can be used to 
select the specific black hole state to generate and includes all possible decay 
modes properly weighted by conserving SU(3) C and U(l) cm . In the following, 
we refer to these different black holes states as quantum black holes, with 
the understanding that we are considering a model which is an extrapolation 
from the classical regime All the QBH code, including PYTHIA 8, is 
written in C++. 



2. How to use the Generator 

The generator consists of four C++ classes, the main class QuantumBlackHole, 
a derived LHAup user class dLHAup, and two utility classes Decay and Random. 
The user normally only needs to interface to the QuantumBlackHole class, 
and only for initialization. 

The package is typically used as two files: the header file qbh . h containing 
class definitions and the file qbh.cc containing the bodies of the member 
functions. 

The user must supply his own main program to initialize the package 
and generate events. The procedure is very similar to writing a standalone 



3 



PYTHIA 8 main program. Above the main program, the following lines must 
appear in this order. 

#include "Pythia.h"; 
using namespace Pythia8; 
#include "qbh.h"; 
using namespace QBH; 



Near the top of the main program, the three objects QuantumBlackHole, 
dLHAup, and Pythia need to be allocated. The following is an example of 
creating them dynamically. 

Pythia* pythia = new PythiaO ; 
dLHAup* lhaPtr = new dLHAup () ; 

QuantumBlackHole* qbh = new QuantumBlackHole (pythia, true) ; 

The second argument to the QuantumBlackHole constructor needs to be true 
for initialization. After initialization, the QuantumBlackHole object may be 
deleted. 

The user is free to further initialize any QBH or Pythia variables before 
the final initialization of the generator. The final initialization of the gener- 
ator is made by passing a pointer to the derived class dLHAup to PYTHIA as 
follows. 



bool status = pythia->init (lhaPtr) ; 



After this, the main program looks like any other PYTHIA main program. 
The user is free to generate events and study these events using all the 
machinery of the PYTHIA classes. 



Les Houches accord event files 13] may be written using the derived class 



dLHAup member functions openLHEF, initLHEF, eventLHEF, and closeLHEF. 

2.1. User Initialization 

The QBH package initializes automatically with a sensible set of parame- 
ters. But typically the user will want to change a few parameters to set up the 
run of interest. The parameters are available as static global variables of the 
QuantumBlackHole class. The first time the QuantumBlackHole constructor 
is call with a boolean argument that is true, a more detailed initialization 
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occurs. Any further instantiations of the class, with or without arguments, 
does no further initialization. However, the user is free to initialize any of 
the parameters using a collection of "set" methods. 

The constructor sets the lower bound on the quantum black hole mass to 
the Planck scale and the upper bound on the mass to three times the Planck 
scale or the proton-proton centre of mass energy, whichever is less. The 
constructor also initializes the LHAPDF parton distribution functions, which 
are the default for QBH, and calculates the branching ratios. Lower bounds 



on the cross section using the calculations of Yoshino and Rychkov [14| may 
be used. In this case, the constructor reads in the appropriate data file. 

Normally the user uses the "set" methods to set the quantum black hole 
state, Planck scale, number of dimensions, etc. For example, the quantum 
black hole state is set using the member functions setQstate, setlstate, 
and the Planck scale is set using setMplanck. Table [1] shows a complete set 
of "set" methods available to the user. 

The method setMaxmass, along with setMinmass, allow the user to set 
a range of black hole masses. This could be useful for restricting the mass 
range for quantum black holes and thus cutting out semiclassical black hole, 
sampling the differential cross section by using narrow mass bins, or for user 
debugging purposes. If the user sets the maximum black hole mass above the 
centre of mass energy, set by method setEcm, the program at initialization 
will reset the maximum black hole mass to be the centre of mass energy. 

There are a corresponding set of methods which return the values of the 
parameters: ecm(), mplanckO, minmassO, maxmassO, qstateO, istateO, 
totdimO, planckdef () , lhaglueO, qscaleO, yrformO, trapO, poissonO, 
RS1(), sm(), chiralO, majoranaO, higgs(), and graviton(). 

It is necessary to define some of the parameter here so their meaningis 
clear to the user. A detailed description of the model can be found in Ref. 
The mass scale used by the Randall-Sundrum black hole (setPlanckdef (1)) 
is defined as 

& = a) 

where m\ is the first KK graviton mass, x\ = 3.83, and c = k/Mp\, where 
k is the ADS curvature and Mpi is the reduced Planck scale. In ADD, the 
fundamental Planck scale in the Dimopoulos-Landsberg, PDG, and Giddings- 
Thomas definitions are, respectively, 
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Table 1: QuantumBlackHole class public member functions for setting parameters. The 
default values are shown as (D=...). 



setQscale (bool) : 



setYRf orm(bool) : 

setTrap(bool): 

setRSl(bool): 



setSM(bool): 
setMajorana(bool) : 



setChiral (bool) : 



setHiggs(bool): 
setGraviton(bool) : 
setTotdim(int) : 

setPlanckdef (int) : 



setLHAglue (int) : 



setQstate(int): 



setlstate(int): 
setMplanck (double) 
setMinmass (double) 
setMaxmass (double) 
setEcm (double) : 



(D=true) definition of QCD scale for PDFs. 

= false: quantum black hole mass, 

= true: inverse gravitational radius. 

(D=false) use Yoshino-Rychkov factors [14] . 

(D=false) use Yoshino-Rychkov trapped surface calculation [14] . 

(D=false) RS or ADD black hole. 

= false: ADD black hole, 

= true: Randall-Sundrum type-1 black hole. 

(D=true) conserve global symmetries. 

(D=false) neutrinos are Majorana or Dirac. 

= false: Dirac, 

= true: Majorana. 

(D=true) handedness of neutrinos. 

= false: neutrinos are left-handed only, 

= true: neutrinos left-handed and right-handed. 

(D=true) include Higgs as allowed particle. 

(D=true) include graviton as allowed particle. 

(D=10) total number of spacetime dimensions. 

= 5 — 11 allowed. 

(D=3) definition of Planck scale [15]. 

= 1: Randall-Sundrum, 

= 2: Dimopoulos-Landsberg definition, 

= 3: PDG definition, 

= else: Giddings-Thomas definition. 

(D=10042) LHA glue code [16]. 

= 10042: LHA cteq611.LHpdf, 

= 19050: LHA cteq5m.LHgrid, 

= 29041: LHA MRST981.LHgrid, 

= else: CTEQ5L (internal PYTHIA 8). 

(D=9) three times electric charge. 

= 4, 3, 2, 1, 0, -1, -2, -3, -4 allowed, 

= 9 used for all partons. 

(D=3) initial state. = 0: q-q, = 1: q-g, = 2: g-g, = 3: all. 
(D=1000.0 GeV) fundamental Planck scale. 
(D=1000.0 GeV) minimum quantum black hole mass. 
(D=14000.0gGeV) maximum quantum black hole mass. 
(D=14000.0 GeV) proton-proton centre of mass energy. 
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where Gd is the D- dimensional Newton constant in the Myers-Perry metric. 

2.2. User Process Run Information 

The PYTHIA initialization call pythia->init (lhaPtr) calls the virtual 
function setlnitO of the derived class dLHAup to initialize the physics 
process in the QuantumBlackHole class. Virtual function setlnitO calls 
QuantumBlackHole : : init () to calculated the total proton-proton cross sec- 
tion, and then QuantumBlackHole: : banner (lhaPrt) to write a banner to 
the output. The banner lists all the parameters, the beam, and quantum 
black hole configuration, as well as, the cross section and estimated statis- 
tical error. Some simple consistency checks are made, and the user may be 
warned or asked to supply different parameters. 

The internal processing of the run information is described next. Al- 
though the user normally does not touch this part of the class, it is informa- 
tive to know how the Les Houches Accord is implemented. The generator 
assumes both beams are protons to use the symmetry of the collision to gain 
some efficiency The LHA weighting strategy is declared by setting the mas- 
ter switch IDWTUP to +1. The generator supplies positive weighted events, 
and PYTHIA returns events with weight 1 by using an accept or reject al- 
gorithm. Only one physics process is declared in the interface to PYTHIA. 
The selection of which quantum black hole to generate is determined by QBH 
according to the user settings. The user may select any of the 14 quantum 
black holes states for generation, or a mixture of all the states weighted by 
their relative cross section is possible (see Table [2]). 

2. 3. User Process Event Information 

The user generates an event by calling pythia->next () which calls the 
virtual function setEventO of the derived class dLHAup to generate one 
event with the QuantumBlackHole class. Virtual function setEventO calls 
QuantumBlackHole: : event () to generate an event. 

The differential cross section is integrated and sampled using importance 
sampling. The integration and importance sampling is very fast since the 
cross section is only a function of two variables and transformations are 
applied to flatten the function that is sampled. 
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Table 2: Qstate and Istate codes for the 14 quantum black hole states. 



State 



Qstate Istate 




Q BH dd 

Q BHl J 9 
QBH x J d 

QBH° qq 

QBH° gg 



4 



3 



2 



2 



























2 





-1 





QBH dg 1/3 


-1 


1 


QBH M 2/3 


-2 





QBH^ g 2/3 


-2 


1 


QBE ~ d \ 3 


-3 







-4 





all 


N/A 
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Sometimes the sampled cross section is vanishingly small. This can occur 
near the edges of the parameters space or if the trapped surface cross section 
is being calculated. In these cases, the partons are both set to up-quarks 
and setEventO still returns true. If false was returned, the cross section 
calculated by PYTHIA would not be the same as that calculated by QBH. 

3. QuantumBlackHole Class Member Functions 

Figured] shows the program flow through the QuantumBlackHole member 
functions. Each public member function is described next. 

bool dLHAup::setInit(void) 
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QBH::setState 



LHAup::initLHEF 



Pythia::next 



Pythia::eventLHEF 
i 



LHAup::setEvent 



Pythia::statistics 



QBH::event 



LHAup::setProcess 

i 



LHAup::setPdf 
1 



QBH::decprd 
I 



QBH::flow 
1 



QBH::deckin 



I 



LHAup::addParticle 



QBH::xpart 
i 



QBH::pdf 
"I 



QBH::xsec 
i 



QBH::sample 



Pythia::init 




LHAup::Setlnit 





LHAup::setBeam 



LHAup::setStrategy 

i 



QBH::init 



LHAup::addProcess 



QBH::banner 



QBH::xpart 
I 

QBH::pdf 
I 

QBH::xsec 



Figure 1: Program flow through the member functions, xpart, pdf , xsec, and sample are 
private member functions. 



Purpose: 
Objects declared 
Functions called: 

Return: 



Les Houches initialization virtual function. 
QuantumBlackHole. 

LHAup::setBeamA, setBeamB, setStrategy, addProcess, 
QuantumBlackHole::init, banner, 
always true. 



setlnit sets the beams characteristics, declares the generation strategy, and 
defines a single physics process. It also calculates the total proton-proton 
cross section and prints a program banner on the standard output. 

bool dL HAup:: set Event (int idProcess) 
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Purpose: 
Objects declared 
Functions called: 



Les Houches event virtual function. 



Return: 



QuantumBlackHole. 

LHAup::setProcess, setPdf, addParticle, 
QuantumBlackHole: :event, decprd, flow, deckin. 
always true. 



setEvent sets the weight for each event, sets some PDF information, and 
adds the beam and generated hard scattering particles to the event record. It 
also adds the quantum black hole state to the event record for documentation 
and process identification purposes. It generates the event by calculating a 
weight for the event and determines random Q, x and beam particle IDs. It 
determines the quantum black hole decay products, connects the colour flow, 
and determines the mothers of the partons. The kinematics of the decay are 
calculated. 

QuantumBlackHole:QuantumBlackHole(Pythia*,bool initialize) 
Purpose: QuantumBlackHole constructor. 

Functions called: Many set methods. 

Return: All user parameters, initializes PDFs, and calculates branching 



ratios. 



The constructor sets information about the quantum black hole and extra 
dimensions. These parameters should be fixed for the entire execution of the 
program. 



void QuantumBlackHole: :init(double& dxsec) 

Purpose: Quantum black differential cross section calculation. 

Functions called: QuantumBlackHole: :xpart, pdf, xsec. 
Return: proton cross section (pb) at random x values. 



init calculates the parton differential cross section at a random x m i n . It 
then determines two random x values and stores the PDFs. It calculates the 
proton-proton differential cross section at these x values. 



void QuantumBlackHole: :banner(dLHAup* inf) 
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Purpose: Output initialization information. 

Functions called: none. 
Return: none. 

banner outputs initialization information in a banner. 

bool QuantumBlackHole::event(double& wxsec,double& Q, double* xx,int* 
idn) 

Purpose: Quantum black hole event generation. 

Functions called: QuantumBlackHole::xpart, pdf, xsec, sample. 
Return: always true. 

event calculates the parton cross section at a random x min value and returns 
the resulting differential cross section and Q value. It then determines two x 
values and uses these to store the PDFs. The proton-proton cross section is 
sampled to determine an event weight and the two partons that caused the 
process are returned. 

void QuantumBlackHole::decprd(int *idn) 

Purpose: Determine decay products identifiers. 

Functions called: ParticleDataTable::chargeType. 
Return: PDG ID codes of decay particles. 

decprd determines the decay products based on predetermined branching 
ratios. 

void QuantumBlackHole::flow(int *idn,int moth [5] [2], int icol[5][2]) 

Purpose: Determine the mothers and colour connections. 

Functions called: none. 

Return: Mother and colour connections. 

flow sets the mothers and colour connections for each final-state particle. 
In most cases, the colour flow is straightforward. Figure |2] shows the set of 
colour flow diagrams for most of the quantum black hole decay processes. 

Figure [3] shows the two diagrams in which the colour flow can not be 
implemented naturally, and we have used a non-standard procedure. The 
technique is as follows. We replace one of the incoming quarks or antiquarks 
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Figure 2: Colour flow diagrams for most the quantum black hole decays. The symbol X 
represents a colourless particle. 

by an antiquark or quark, respectively, of the same flavour. The flavour is 
changed just before the particle is added to the event record. Thus the par- 
ticle kinematics and weight of the event are unaffected. Next, an additional 
outgoing quark or antiquark is added so that electric charge is conserved. 
The colour connections can then be made between the coloured particles. In 
this way we can conserve electric charge and allow colour to flow. The addi- 
tional outgoing particle belongs to the first quark generation and is given zero 
momentum-energy. The extra zero momentum-energy particle adds little to 
the event. Typically a colour string will form between it and a beam- remnant 
diquark. String fragmentation will give a few extra baryons and mesons. 
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Because of the vanishing momentum-energy, the diquark will dominate the 
kinematics and typically the extra particles will travel down the beam pip e 
undetected. Such an approach has been used in previous generators 17|, |l8 . 




Figure 3: Diagrams for which a nonstandard colour flow technique has been used. The 
symbol X represents a colourless particle. 

PYTHIA 8 provides a junction class to deal with baryon violating decays. 
Currently it handles only the case of one initial state particle. The facility is 
not readily available to the users via the Les Houches Accord interface. 

void QuantumBlackHole::deckin(int *idn,dLHAup* lha,double *pl, double *p2) 
Purpose: Calculate kinematics of decay. 

Functions called: Decay: :twobody. 

Return: four-vectors of black hole and final-state particles. 

deckin calculates kinematics of the decay. 

The following four member functions are private. 

void QuantumBlackHole::xpart(double& xmin,double& Q,double& sighat) 
Purpose: Calculate parton-parton cross section (pb). 

Functions called: none. 
Return: xmin, Q, sighat. 

xpart generates a random x min and uses this to calculate the parton-parton 
differential cross section. It also returns Q which is only available after the 
gravitational radius and mass are calculated. 

void QuantumBlackHole::pdf (const double* xx, const double Q) 

Purpose: Get parton distributions via PYTHIA PDF class. 

Functions called: PDF::xf. 

Return: fills disf[2][13] array. 
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pdf gets the parton distribution functions via the PYTHIA PDF class and 
stores them as a private data member of the class. The PDFs are available 
to this function via a stored private pointer in the class. 

void QuantumBlackHole::xsec(const double sighat,double& dxsec) 

Purpose: Calculate the proton-proton cross section (pb). 

Functions called: disf [2] [13] . 

Return: Differential cross section in pb 

xsec uses the parton cross section and an array of PDFs to calculate the 
proton-proton cross section. 

void QuantumBlackHole::sample(const double sighat, const double dxsec, int* 
idn) 

Purpose: Sample cross section to select beam partons. 

Functions called: disf [2] [13] . 

Return: weight, beam parton IDs. 

sample samples the cross section to determine the beam partons and returns 
a weight. 

3.1. Utility Classes Random and Decay 

The totally static Random class takes care of all the random distributions 
and random selection utilities, such as selecting a random charged lepton, 
for example. The class is driven by the uniform random number generator 
used by PYTHIA. Thus PYTHIA may be used to initialize and manage the 
seeds of the Random class. 

The Decay class handles the two-particle phase space decay. It has one 
public member to perform the decays. A set of private member functions 
perform all the required Lorentz transformations. 

4. Installation and Availability 

The QBH generator is written entirely in C++. It interfaces only to 
PYTHIA 8. If a lower-bound on the cross section is being calculated by using 
the trapped surface option, then the appropriate file is read in. These small 
text files must exist in the same directory as as the executable image. If the 
LHAPDF parton densities are required, the LHAPDF library must be linked. 
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Figure 4: Total cross section versus number of loop iterations used in the cross section 
calculation. 



The generator has been tested with PYTHIA version 8.150 and LHAPDF 
version 5.8.0. Version 1.02 and earlier of the generator should be used with 
version 8.130 of PYTHIA. Both Linux and Darwin-gcc4 architectures with 
the slc4Ja32_gcc34 platform have been tested. 



The source code can be obtain from the HepForge site http: / / projects.hepforge.org/ qbh/ 



5. Performance 

The performance of the event generator has been measured. Figure H] 
shows the cross section versus number of loop iterations use to calculate 
the cross section for three different runs at each value of the number loop 
iterations. The spread in the three cross section values at each number of loop 
iterations is consistent with the cross section statistical error bars. The cross 
section calculation is stable and converges nicely. For a statistical accuracy 
of 0.2%, 10 5 iterations are required, while for an accuracy of 0.07%, 10 6 
iterations are required. The totally inelastic cross sections have been used 
for the performance studies. One might expect the trapped surface cross 
sections to take a little longer. 
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As a further check of the cross section calculation, we compared the cross 
section calculated by QBH using 10 5 loop iterations with that calculated by 
PYTHIA for 10 5 generated events. The percentage difference was less than 
0.5%. 

Figure [5] shows the processor time versus number of loop iterations in the 
cross section calculation. To obtain a 0.2% accuracy requires 2 s, while for a 
0.07% accuracy 14 s is needed. However the processor used was arbitrary so 
the time axis should be viewed as a relative time. By default, the generator 
used 10 5 loop iterations to calculate the cross section. 




j_d i i i i i i i i 

10 4 1 5 ,10 6 . 10 7 1 0' 

Loop Iterations 



Figure 5: Processor time versus number of loop iterations in the total cross section calcu- 
lation. The timing resolutions is 1 s and the processor time axis should be considered as 
a relative timing. 

Figure [6] shows the processor time versus number of generated events. To 
decouple the generator from PYTHIA 8, the "hard process" times have frag- 
mentation, hadronization, particle decays, initial- and final-state radiation, 
and multiple interactions switched off. About 6 x 10 3 QBH events can be 
generated per second. If fragmentation, hadronization, and decays are turned 
on, about 1 x 10 3 events can be generated per second. Turning on initial- 
and final-state radiation allows 30 events per second to be generated. Finally, 
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Figure 6: Processor time versus number of events generated. The timing resolutions is 1 s. 
The processing tasks are cumulative. For example, multiple interactions also includes all 
the other processing. 



turning on multiple interactions (everything) allows 8 events per second to 
be generated. 
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